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Abstract 

We study the production amplitude for the reaction NN — > NNir up to next-to-leading order in 
chiral perturbation theory using a counting scheme that takes into account the large scale introduced 
by the initial momentum. In particular we investigate a subtlety that arises once the leading loop 
contributions are convoluted with the NN wavefunctions as demanded by the non-perturbative 
nature of the NN interaction. We show how to properly identify the irreducible contribution of 
loop diagrams in such type of reaction. The net effect of the inclusion of all next-to-leading order 
loops is to enhance the leading rescattering amplitude by a factor of 4/3, bringing its contribution 
to the cross section for pp — > dir + close to the experimental value. 



1 Introduction 



The highly accurate data for pion production in nucleon-nucleon collisions close to the production 
threshold are a challenge for theoreticians. When the first close to threshold data for the total cross 
section of the reaction pp — ► ppir appeared in 1990, existing models fell short by a factor of 5-10 
^H]. Many different mechanisms were proposed to cure this discrepancy: heavy meson exchanges 
[3] , (off-shell) pion rescattering jlj , excitations of baryon resonances [H] , and pion emission from 
exchanged mesons [7j- The total cross sections for the reactions pp — ► pmr + and pp — ► dir + on the 
other hand could always be described within a factor of 2 — the amplitude is dominated by the 
isovector rescattering contribution 1 . For a recent review see Ref. (Hj. 

To resolve the situation various groups started to investigate NN — ► NNtt using chiral perturbation 
theory (ChPT). As an effective field theory (EFT) it is to be free of any ambiguities and people 
expected that now the relevant physics of NN — ► NNtt could be identified. As a big surprise to 
many, however, it turned out that, when naively using the original power counting by Weinberg jX, 
the discrepancy between theory and data became even larger at next-to-leading order (NLO) for 
pp — ► pp7T° PHI as well as for pp — > dir + . At the same time it was already realized that a modified 
power counting is nescessary to properly take care of the large momentum transfer characteristic for 
pion production in NN collisions, however, when applied in its original formulation the results were 
basically the same for neutral ^2] as well as charged pions ^Sj. Even worse, the corrections at one- 
loop order (next-to-next-to-leading order (N 2 LO) in the standard counting) turned out to be even 
larger than the NLO corrections, indicating a divergence of the chiral expansion [141 115j . 

Recently there were two developments: one that focused on some technical aspects related to the 
evaluation of the matrix elements |16l I17j and another regarding the power counting for the large 
momentum transfer reactions. Formal inconsistencies of the naive power counting using the heavy 
baryon scheme were pointed out in Ref. [T^]. I n addition the ideas formulated in Refs. ^3 
were further developed and improved — it was especially recognized how to properly estimate loop 
contributions. This scheme was implemented in Refs. |19| I20| — the essential features are described in 
detail below. The basic conclusion was that an ordering scheme exists for the reactions NN — ► NNtt 
that can lead to a convergent series. However, so far full calculations (including the distortions due 
to the N N interactions) within this scheme exist only for the production of p-wave pions ^S] • 

In Ref. [201 it was demonstrated by explicit evaluation of the leading loop contributions (shown in 
figure mb)-(d)) how the presence of the large momentum scale influences loops. The central findings 
of that work were that it is possible to define an ordering scheme for NN — > NNtt, but some loops are 
to be promoted to significantly lower orders compared to what is expected from Weinberg's original 
counting. Ref. [201 left two questions unanswered that we will address in this note: 

• For the reaction channel pp — ► ppir the sum of all leading loops canceled; the origin of this 
cancellation could not be identified; 

• For the channel pp — * dir + the sum of diagrams of Fig. H^b)-(d) gave a finite answer. However, as 
pointed out recently in Ref. [2J, the corresponding amplitudes grow linearly with increasing final 
AW-relative momentum. This behavior leads to a large sensitivity to the final NN wavefunction, 
once the convolution of those with the transition operators is evaluated as demanded by the non- 
perturbative nature of the NN interaction. The solution to this problem proposed in Ref. [2^ 
is to include a new counter term at leading order to absorb this unphysical behavior. However, 
chiral symmetry does not allow for such a structure (see Appendix for details). 

As we will show in this paper, the solution to both questions are related and at the same time shed 
some light on the concept of reducibility in pion reactions on few-nucleon systems. 
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Figure 1: Leading loop diagrams for NN — > NNir. Here dashed lines denote pions and solid lines 
denote nucleons. 

We further demonstrate that the net effect of the inclusion of the NLO loops, shown in Fig. ^ is to 
enhance the leading rescattering amplitude by a factor of 4/3, bringing its contribution to the cross 
section for pp — > c?7r + close to the experimental value. 

2 Power counting and the concept of reducibility 

Already the existence of nuclei shows that perturbation theory is insufficient to properly describe 
two-nucleon systems: only an infinite sum of diagrams can produce a pole in the S'-matrix. To bring 
this observation in line with power counting, Weinberg proposed to classify all possible diagrams 
according to the concept of reducibility |221 1231 124j : those diagrams that have a two nucleon cut are 
called reducible. Those which do not are called irreducible. The latter make up the potential that 
is to be constructed according to the rules of ChPT. The former are then generated by solving the 
Schrodinger equation, using the mentioned potential as kernel. This scheme acknowledges that the 
two nucleon cut contributions are enhanced compared to the irreducible parts. 

It was also Weinberg who gave a recipe how to calculate processes on few nucleon systems with 
external probes [S]: here the transition operators are to be calculated using ChPT. Then those transi- 
tion operators must be convoluted with the appropriate AW wave functions — in full analogy to the 
so-called distorted wave Born approximation traditionally used in phenomenological calculations pQ. 

Therefore it is necessary to disentangle those diagrams that are part of the wave function from those 
that are part of the transition operator. In complete analogy to AW scattering described above, the 
former are called reducible and the latter irreducible. Also here the distinction stems from whether 
or not the diagram shows a two nucleon cut. Thus, in accordance to this rule, the one loop diagrams 
shown in Fig. ^b)-(d) are irreducible, whereas diagrams (a) seem to be reducible. However, it will be 
the central finding of this work that diagrams (a) contain a genuine irreducible piece due to the energy 
dependence of the leading AW7T7r-vertex, the so-called Weinberg-Tomozawa term (WT). Specifically, 
the energy-dependent part of the WT vertex cancels one of the intermediate nucleon propagators, 
resulting in the irreducible part of diagrams (a). 

As mentioned in the introduction, the power counting needs to be modified in order to be applicable 
for AW — > NNir. The reason for this necessity is the magnitude of the nucleon center-of-mass 
momentum p required to produce a pion at rest in AW collisions. It is given by 



with M = 939 MeV and m n = 139.6 MeV denoting the nucleon and pion mass, respectively. Eq. (|T|) 
exhibits the important feature of the reaction AW — > NNtt, namely the large momentum mismatch 
between the initial and the final nucleon-nucleon state. This leads to a large invariant (squared) 



|p| = yW(M + m^/A) 



(1) 
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Figure 2: Tree level diagrams that contribute at leading ((a) and (b)) and next-to-leading order (c) 
to AW — > NNtt. The double line denotes a A-isobar. Note, in diagrams (b) and (c) — for illustration 
— with the one-pion exchange only one part of the AW and AW — > NA potential is shown. 



momentum transfer t = —Mm n between in- and outgoing nucleons. The appearance of the large 
momentum scale yJMm^ in pion production demands a change in the chiral power counting rules, as 
pointed out already in Ref. ^3 El]. In addition, it seems compulsory to include the Delta-isobar as 
an explicit degree of freedom, since the delta-nucleon mass difference A = 293 MeV is comparable to 
the external momentum p ~ \/Mm^ = 362 MeV. The hierarchy of scales 



suggested by this feature is in line with findings within meson exchange models where the Delta-isobar 
gives significant contributions even close to the threshold |25l 126) ff 1 . The natural expansion parameter 
therefore is 



As a result at leading order only tree level diagrams contribute to the transition operator (diagram (a) 
and (b) of Fig. |2J). Already at next-to-leading order — in addition to the first diagram that involves 
a A-isobar (diagram (c) of Fig. [2j) — the first loops appear (see Fig. As a consequence of the two 
scales p and given in Eq. (J2J) there exists a dimensionless parameter that is of order namely 
m^/p, that can appear as the argument of non-analytic functions as a result of the evaluation of loop 
integrals. Thus, each loop now contributes not only to a single order, but to all orders higher than 
the one where it starts to contribute || . In this work we only consider the leading parts of the loops 
in Fig. [J that start to contribute at NLO. 

At threshold only two amplitudes are allowed to contribute to the reaction AW — > NNtt, namely 
An an d Aio, where we used the notation A^Tr to label the total isospin of the initial (Tj) and final 
(Tf) AW-pair. The third amplitude allowed by the standard selections rules — Aqi — has to have 
at least one p-wave in one of the final subsystems. To the reactions pp — > ppn° and pn — ► ppir~ only 
An and to AW — ► dir only Aio contribute at threshold, whereas both An and Aio contribute to the 
reaction pp — ► pnir + . 

The only transitions that are allowed to contribute near threshold are 3 Po — * 1 Sqs for An and 
3 Pi — ► 3 Sis for Aio, where small letters denote the pion angular momentum with respect to the AW 
system and the AW partial waves are labeled with the standard notation 2S+1 Lj. Those lead to the 
following amplitude structures |Hj 



* For the channel pp — > ppir strong support for an important role played by the A-isobar was given by a partial wave 
analysis |27|. 



m n <C p ~ A < M 



(2) 




(3) 




(4) 



with %' = (xl^yiXs^J /V2 and S = (xT Vy&Xz) / and = (xi^yixlV)/^ where the Xl,2 
(X3,4) denote spinors for the incoming (outgoing) nucleons. For a deuteron in the final state we need 
to use 

M =A w ((Sxp). eA , (5) 



where now * denotes the deuteron polarization vector and A\o the convolution of the production 
operator with the deuteron wave function — although the transition operators for Aiq and Aiq are the 
same, the matrix elements in general have different units. We come back to the evaluation of the 
deuteron cross section in section |1J The amplitude in Eq. @ is normalized such that for the total 
cross section for pp — * pnir + we have 

a tot (PP^ P nn-) = J% [ (\A U \ 2 + 2\A W \ 2 ) , (6) 



lGvr 3 ^ , 

where q' denotes the cms momentum of the outgoing it, p' the relative momentum of the outgoing pn 
pair, s the invariant energy of the system, and u> q = \/ q 2 + m 2 . Analogously, one finds for the total 
cross section for pp — ► dix^ 

+ 2M d M 2 , ~ 2 
ototipp d7T ) = (IP A io , (7) 

TTS 

where M d denotes the deuteron mass. 

If we neglect all NN distortions, we get for the leading rescattering contribution (Fig. Eta)) at 
threshold [20] 

^-(t)^j(f)^-|i( 1+0 (x,), 

where we used that the irN — * irN amplitude to leading order contains not only the standard WT 
term that scales as the sum of the incoming and outgoing pion energies, here equal to m n and 
respectively, but also its recoil correction equal to p 2 /2M = m n /2. The relevant terms of the un- 
derlying Lagrangian density are given in the Appendix. Note that with a value of 2m- K /4:f 2 the WT 
vertex including the recoil correction as it appears in the iViV — * NNir amplitude takes exactly the 
value it has for elastic ttN scattering at threshold. 



3 Evaluation of loops 

In the reaction AW — > NNw the energies of the initial nucleons are of order p 2 /2M ~ and the 
momenta are of order p. In irreducible loops, on the other hand, both energies and momenta are of 
order p (see Appendix E of Ref. |H] for details). As a consequence in the evaluation of diagrams (b),(c), 
and (d) of Fig. ^the nucleon recoil terms (i.e. the nucleon kinetic energies) can be neglected in the 
nucleon propagators, since they scale as m n . At first glance this seems to be at variance with the recent 
finding that three body NNir cuts that originate from the nucleon recoils play an essential role in pion 
reactions on few nucleon systems |28U29j . However, the reactions studied in these references had very 
different kinematics, as they were small momentum transfer reactions, where in typical kinematics the 
ttNN state was near on-shell. Here, on the other hand, we are faced with a large momentum transfer 
reaction: for typical kinematics a ttNN intermediate state is far off-shell. In addition, here there is an 
additional kinematical suppression for the irNN cuts: at the cut the typical pion momentum, which 
set the scale for the typical loop momentum, is that of the external pion of order of at most m^. As 
a consequence the ttNN cuts do not contribute before N 5 LO to the reaction iV./V — ► NNtt |S]. 
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By comparison to the full results of Ref. Jl] , in Ref . |2U| it was shown that it is allowed to expand 
the integrand of the loop integrals before evaluation in powers of ^Jm^/M . As a consequence it was 
possible to express the leading contribution of all loops corresponding to diagrams (b)-(d) of Fig. ^ 
in terms of a single integral 

, ~2s = l [ J^l 1 (q\ 

0{P0,P j ij {2nY{l -ieW-ml){{l+pY-mD ' { ' 

One finds I (Po,p 2 ) = h{p 2 ) (1 + O (x)) = 1/(16^^2) (1 + (x)). The assumption of threshold 
kinematics (all outgoing momenta vanish) simplifies the operator structure significantly and we can 
write — neglecting for the moment the distortions from the NN interaction — in order [2U] 

A\ b + lc+ld = _^(_2 + 3 + 0) p 2 / (p 2 )-^^ 



16/| v 256/ s 

A U +1C+W = ^(-2 + 3-1) p 2 / (p 2 ) = 0. (10) 

Note, here and in what follows we write equalities although we dropped terms of higher order in \- 
As mentioned in the introduction, the sum of the NLO loops vanishes in case of An. We will give 
an explanation for this cancellation below. Let us now concentrate on A\q. In order to compare 
the result of Eq. (|10|) to data, the transition operators need to be convoluted with appropriate NN 
wavefunctions. The convolution integrals that arise necessarily involve non-vanishing p' , denoting the 
outgoing NN relative momenta, even if we still work at threshold. However, the structure of the loop 
integral Iq is such that — to leading order — only p — p' appears in the integrals and thus one can 
directly generalize the expressions of Eqs. (jTU|) . As was argued in Ref. [25 > this implies that for large 
p' the contributions from the loops grows linearly with \p'\. When it then comes to the convolution of 
those operators with NN wave function this linear growth of the transition operators leads to a large 
sensitivity to the deuteron wave functions. However, there should be no sensitivity to the particular 
wave functions used, for off-shell quantities are not observable |30| I31j. On the other hand, the chiral 
Lagrangian does not allow for a counter term to compensate this linear growth — see Appendix for 
details. The solution given in Ref. (2^, namely the inclusion of a counter term at leading order, is 
therefore not consistent with the effective field theory used. However, as we will show, the loops with 
the unwanted behavior will be canceled exactly by an irreducible piece of diagrams (a) of Fig. We 
proceed as follows: we first show this cancellation to one loop order without distortions. Then we 
generalize the result to the inclusion of the full NN wave functions. 

We still assume threshold kinematics and now turn to the evaluation of diagrams (a) of Fig. ^ In 
doing so one first has to realize that in contrast to the irreducible diagrams discussed in the beginning 
of this section, energies in the diagrams with a two-nucleon cut are of the order of the external energies 
{Iq ~ p 2 /2M ~ m n ). Therefore, there is a priori no reason to neglect the nucleon recoils that are of 
the order of m n . We thus get for the full expression for the first diagram of Fig. ^a), up to higher 
orders, 

A lal_^9\ f dAl [l + m w -(2p + l)-l/(2M)} 



10 " 32/1 J { 2,y (lo _ ^ _ + . e)Ho + ^ _ + ie) 

(i -(Up)) 



(p-mDdl+pY-ml) ' 
A\f=0, (11) 

where we included the recoil correction to both the WT term in the numerator as well as the nucleon 

llal 

hi 



energies in the denominator, in line with the discussion above. The vanishing A\^ reproduces the 



6 



well known result that the WT interaction does not contribute to the leading rescattering diagram in 
pp — > ppir . 

In order to proceed we rewrite the first term in the numerator of the above integral as 



l + m n 



(2p + I) ■ I 
2M 



2M 



2m„ 



where we used that at threshold p 2 /M = m w . The first term now exactly cancels the first nucleon 
propagator and we are left with an expression that no longer has a two nucleon cut — it is irreducible. 
In this irreducible piece we can neglect the recoil corrections in the remaining nucleon propagator - 
c.f. the discussion at the beginning of this section — and get 



ilol 
l 10 



+ 



■ 3ffj f dH J (1-(1+P)) 

''32/5 J (27r) 4 | ( _ /o + ie)(Z 2_ m 2 )((Z+p) 2_ m 2) 

2m^ (f. (f+p)) 



(lo 



2 



(I+P) 2 
2M 



(i+p) 2 

2M 



(I 2 — m 2 ){(l + p) 2 — m 2 ) 



(12) 



Up to higher orders the first term gives 



.lai(irr) 
^10 



9A ~2 



4y 16/5 



3 g\\P\ 

4 256/5 ' 



(13) 



where the label (irr) indicates that this is only the irreducible piece of the diagram. Analogous 
considerations for the second diagram of diagrams (a) of Fig. ^gi ve 



\la2 

L io 



dH 



32/|7 (2VT)4 |(-/ + i e )(/2_ m 2 )((/+p) 2_ m 2 ) 

2m^ (l-(l+p)) 



2 2M 

The leading piece of the first term gives 



A 



ia2(irr) 
10 



9 A „,2 



4 J 16/5 



P hip' 



1 9a\p\ 
4 256/5 • 



Thus we get 



A 



iai(irr)+ia2(irr)+ife+ic+id 
10 



A iai(irr)+ia2(irr)+ifc+ic+id 



16/| 
16/1 



2 + 3 + ) p 2 I (p 



( 0+0-2 + 3-1) p 2 I (p 2 ) = 



(14) 



(15) 



(16) 



where we repeat the results for An from above for comparison. Thus, in both channels that contribute 
at the production threshold the sum of all irreducible loops that appear at NLO cancels. On the other 
hand the remaining pieces in the expressions for A\q exactly agree to the convolution of the leading 
rescattering contribution with the one pion exchange, however, with the NNirir WT vertex put on- 
shell. Thus the WT vertex takes the value 2m 7T /(4f 2 ) — c.f. Eq. (jHJ). The two-nucleon propagators 
in these integrals have a unitarity cut and it is this cut contribution that should dominate the integral 
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a) b) c) 

Figure 3: Illustration of the kind of topologies that corresponds to the irreducible structures (denoted 
by the filled box on the propagator that gets canceled by the energy dependence of the irN — > ttN 
vertex) that emerge from the convolution of the energy dependent rescattering term with various 
contributions to the AW potential. 



- in line with Weinberg's original classification as reducible and irreducible. In other words, these 
pieces are indeed dominated by the reducible piece* 2 . 

Next we show that for all ingredients of the AW potential but the one-pion exchange the choice of 
an on-shell NNittt vertex is of sufficient accuracy. To this end we remind the reader that the integral 
corresponding to the irreducible pieces of this first diagram of Figure ^a) had a structure like diagram 
(b) of that figure. This is illustrated in part (a) of Fig. |21 At leading order in the AW potential 
there is besides the one pion exchange also a contact interaction. The convolution of the rescattering 
diagram with the WT vertex (diagram (a) in Fig. [2J) with this part of the AW potential we can again 
decompose into a reducible piece with the ttN — > irN vertex on-shell and an irreducible piece that 
takes a structure of an integral with one nucleon less (see part (b) of Fig. |3J); this diagram, however, 
does not contribute below N 4 LO and is therefore irrelevant to the order we are working. Thus, all 
that is to be kept is the convolution of the on-shell rescattering contribution with the contact AW 
interaction — in line with the findings of Ref. [T^]. At NLO in the AW potential there are pion loops. 
Then the irreducible piece of the convolution of the leading rescattering contribution with this piece 
results in a two loop diagram (for a particular example see part (c) of Fig. |3J) that does not contribute 
up to N 3 LO. Thus, to the order we are working we can safely put the WT vertex on-shell for the 
convolution of any piece of the AW potential with the leading rescattering contribution. 

What remains to be shown is that the cancellation of Eq. (|16|) survives (to the given order) the 
convolution with the full wave functions. This generalization is straightforward. The corresponding 
diagrams are shown in Fig. |1] for the inclusion of the final state interaction that we want to discuss 
in detail. The argument in case of the initial state interaction is completely analogous and will not 
be given. Note that only the irreducible parts of the diagrams (a) are to be included — the reducible 
pieces get absorbed into the wave functions. Let k denote the integration variable of the convolution 
integral that we chose equal to the four momentum of one of the nucleons. As argued in the previous 
paragraph, the integral will be dominated by energies close to the corresponding on-shell energies. 
This sets the scale for /c-especially we can safely assume ko <C p. This is all that is needed to neglect 

in the nucleon propagator of the pion loop integrals. On the other hand in these loops k only enters 
as p — k. Thus, the terms that enter in the convolution integrals with the final state interaction to 
the order we are working are simply given by replacing p by p — k in Eqs. (|10|) . (|13|) . and (|15|) . This 
will obviously not change the relative strength of the individual diagrams — the cancellation survives 

* 2 In addition, the relative strength as well as sign in these terms equal to —3 for the pion exchange in the final Tf = 
channel compared to +1 for the pion exchange in the initial Ti = 1 channel equal to the expectation values of the isospin 
parts of the one pion exchange in the relevant isospin channel, since (TT3\fi ■ T2ITT3) = 2T(T + 1) — 3. 
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a) b) c) d) 

Figure 4: Leading loop diagrams for NN — > NNir, convoluted with the NN T-matrix, denoted by 
the ellipse. Here dashed lines denote pions and solid lines denote nucleons. The filled box on the 
nucleon propagators of the diagrams (a) indicate that only the irreducible piece is to be taken. The 
reducible part gets absorbed into the wave function. 



the convolution with the wave functions. We therefore conclude that up to next -to-leading order all 
irreducible pion loops in the transition operator cancel with the only effect that the WT vertex in the 
rescattering diagram is to be put on-shell. 

4 Results 

Since the NLO diagrams discussed in this work contribute only to A±o, we will now discuss their 
impact on the reaction pp — > dir + that is fully determined by that amplitude. The cross section data 
for this reaction near threshold is traditionally parameterized as 

a = arj + (3rf , (17) 

where Coulomb effects were neglected. To the reaction pn — > dir both a and (3 contribute with only 
half their strength [221 • Here r\ denotes the outgoing pion momentum in units of its mass. The first 
term gives the s-wave strength, whereas the second one denotes the p-wave contribution (as well as 
some possible energy dependencies of the s-wave |2*T]). 

Before comparison with experiment is possible, the transition operators are to be convoluted with 
appropriate AW wave functions. Here we use basically the same formalism as described in Ref. 
29 * 3 and thus we do not give any formulas in detail. For the NN distortions we use the CD-Bonn 
potential [HZj and use the same parameters as in Ref. [22]. To calculate the leading order (LO) 
rescattering process (diagram (a) of Fig. |2J) we use the standard expression for the WT term in 
threshold kinematics — thus we put 3/2711,^ at the vertex In addition we also evaluate the 

direct contribution (diagram (b) of Fig. |2J). Since we derived the transition operator in threshold 
kinematics we can only calculate a. We obtain 

a LO = 131 fib . (18) 

This number is dominated by the rescattering contribution; switching off the direct term alpha is 
lowered by 30 % to a WT ' LO = 101 pb. This value is consistent with those given in Ref. [3S]. Note 
that the direct term is known to be quite model dependent [HE] and is small because of a cancellation 

* 3 For more details see Appendix of Ref. 36 . 
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Ql 1 1 1 1 1 1 

0.05 0.1 0.15 

Figure 5: Comparison of our results to experimental data for pp — > dir + . The dashed and dashed- 
dotted curves show the LO results, where for the latter the direct contribution was omitted. The solid 
and dotted line show the results at NLO, where in the latter the direct term was omitted. The data 
is from Refs. [32^ (open circles), [33] (filled circles) and [31] (filled squares). 



of individually sizable terms. Clearly such a cancellation can not be captured by the counting scheme. 
Still, we point out that in an EFT scheme as used here all terms at a given order have to be retained. 

As outlined above, in order to include the NLO contributions all we need to do is to replace the 
(3/2)77% in the WT vertex by 2m 7r — or, stated differently, to scale the given results for a T ' LO by 
a factor (4/3) 2 . Thus, for the rescattering piece we get at NLO a WT ' NLO = 182 fjb, whereas the full 
result including the direct term is 

a NLO = 220 ^b . (19) 

We checked that using different models for the NN distortions changed the rescattering contribution 
by about 10 %, in line with previous findings [35j . 

In Fig. we compare the results of our calculation to the experimental data* 4 . One clearly sees 
that going from LO to NLO improves the description of the data. Note also that even a change by 
about a factor of 2 in the cross section is in line with what is expected from the counting: after all 
the expansion parameter for the amplitude is x ~ 1/2- 

For a complete NLO calculation what needs to be considered in addition are diagrams of the kind 
depicted in Fig. [2] (c), but with the full NN — > iVA transition amplitude. The evaluation of this 
diagram needs as input the NN — > NA transition potential consistent with chiral symmetry. The 
corresponding parameters are to be fixed from a fit to NN data. Such a fit has not yet been performed 
and thus we postpone the evaluation of diagram (c) to a later work. Note, however, that model 
calculations show that the contribution of the Delta-isobar is not more than 10 % in the amplitude 

m 

* 4 Note that the data from pn — > dn |32| are considerably lower than those for pp — > div + |33l 1341 . It appears unclear 
whether this discrepancy is due to systematic uncertainties (not shown in the figure) or due to other sources |38| . 
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5 Summary and conclusions 



We have shown that the proper set of diagrams that contributes to the transition operator for the 
reaction NN — ► NNir at NLO in chiral perturbation theory is given by the diagrams of Fig. [21 however, 
with the NNirir vertex in diagram (a) put on-shell. To get to realistic results these operators are to be 
convoluted with proper iViV wave functions. The irreducible chiral loops that arise at this order exactly 
cancel those terms that arise from the off-shell parts of the WT vertex. This cancellation is required 
for formal consistency of the whole scheme, since the mentioned diagrams show a linear growth with 
respect to the outgoing NN momentum. Such a growth would have led to a large sensitivity to the 
NN wave function, when the convolution with the final state interaction is calculated. This, however, 
would have been in conflict with general arguments. 

This at the same time also explains, why the sum of all loops has to vanish at next-to-leading 
order for the reaction pp — ► ppir : in this channel there is no leading rescattering contribution. Thus, 
there is also nothing that could cancel the linear divergence discussed above. The consistency of the 
formalism therefore demands the sum of all loops to vanish. 

As a result of our findings we can conjecture a general recipe on how to deal with pion reactions in 
a nuclear environment in the presence of time derivatives in vertices: one has to calculate all diagrams 
up to a given order, including those that are formally reducible. Then the energy dependence in the 
vertices is used to cancel one of the nucleon propagators. This produces an irreducible piece that is to 
be part of the transition operator as well as a reducible piece, where, however, the energy dependence 
of the vertices is replaced by the corresponding on-shell value* 5 . 

This new rule has significant impact on the role of isoscalar rescattering in NN — > NNir (diagram 
(a) of Fig. |2 but with the leading isoscalar interaction used for irN — * irN). Empirically the isoscalar 
irN scattering length is known to be very small. Theoretically it turned out that this smallness is 
a consequence of an efficient cancellation amongst individually large terms |4U| . Due to the energy 
dependence of the irN —* irN operators, however, when evaluated in the kinematics relevant for pion 
production in NN collision this cancellation is much less efficient leading to a significant contribution 
from isoscalar rescattering |12M1U| . If, on the other hand, the above rule is used, isoscalar rescattering 
enters with the strength of the very small isoscalar irN scattering length and thus would give a 
negligible contribution. 

We have demonstrated that the net effect of the inclusion of the NLO loops, shown in Fig. is to 
enhance the leading rescattering amplitude by a factor of 4/3, bringing its contribution to the cross 
section for pp — ► dn + close to the experimental value. 

The next steps will be to evaluate the NN — > NNir amplitudes to N 2 LO for both s- and p-wave 
pions for all possible amplitudes. At this order two counterterms enter for the pion s-waves and one 
for the pion p-waves both accompanied by S-wave nucleons in the final state. To this order p-wave 
pions together with P-wave nucleons are parameter free predictions. On the other hand there are in 
total more than 40 observables* 6 measured for the reaction channels pp — > ppir [41j , pp — > pmr + [12] , 
pp — > dir + [43j and pn — ► ppir~ |44j . Up to now only one phenomenological calculation was compared 
to this large amount of data |45[ [Hj and it was found that all charged channels are well described, 
whereas there are significant discrepancies for the neutral pions. It will therefore be of strong interest 
to see if the new structures that emerge from the chiral Lagrangian are able to cure these discrepancies. 

Once the described channels are analyzed within ChPT one should move ahead to consistently 
investigate the isospin violating observables measured recently, namely the forward-backward asym- 
metry in pn — > dir° 46 as well as the total cross section for dd — ► cot |47j . First steps in this direction 

* 5 Note that in the presence of quadratic time derivatives at individual vertices or the simultaneous appearance of 
several time derivatives in one diagram, additional vertices are to be included |39|. 
#6 This large number is achieved by fully exploiting the 5 dimensional phase space |41|. 
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were taken in Refs. (IHIHU] . 
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A Lagrange densities and vertices 

The starting point is an appropriate Lagrangian density, constructed to be consistent with the sym- 
metries of the underlying more fundamental theory (in this case QCD) and ordered according to a 
particular counting scheme. Omitting terms that do not contribute to the order we will be considering 
here, we therefore have for the relevant terms of the leading and next -to-leading order Lagrangian in 
sigma gauge for the pion field [23] 

+N^[id --^T ■ (tt x ir)}N+^-NW ■ a ■ (v7r+-^7r(7r • Vtt) ) N 

+i^—[N^V 2 N+^ a V 2 ^a] + J ,J iNW ■ (tt x Vtt) • VN+h.c.) 
2m n b%/; 

!]A -[iN^r ■ -ho ■ VN+h.c] ^— [iN*T ■ ivS ■ V^A+h.c] 



4m N f n 2mjv/7r 
-^-N\t . & . V7T ) NN 1 N _ ^ eijkeabcdi7 , aN \ a3TbN N \ akTcN + ... ? ( 20 ) 

where f n denotes the pion decay constant in the chiral limit, is the axial- vector coupling of the 
nucleon. hj± is the ANir coupling, and S and T are the transition spin and isospin matrices, normalized 
such that 

1 1 

j = 3(2% - itijkVk) , T-rj = ~(2<% - ie ijk T k ) . (21) 

The dots symbolize that what is shown are only those terms that are relevant for the calculations 
presented. As demanded by the heavy baryon formalism, the baryon fields ./V and are the velocity- 
projected pieces of the corresponding relativistic fields; e.g. N = 1/2(1+ where = (1,0,0,0) 
denotes the nucleon 4-velocity. The corresponding vertex functions can be read off directly from 
Appendix A of Ref. jSHj- The relevant vertices for the NN interaction are discussed in Ref. |51j . 

The last two terms in Eq. (|20|) show the leading four-nucleon-pion counter terms. The Pauli 
principle for the NN system only allows them to contribute in one fixed linear combination. The 
corresponding operator contributes to pion p-waves in pp — * pnir + as well as to the leading three- 
body force |191 152j . To NN — ► NNtt at threshold this operator can only contribute through loops. 
The leading four-A r -7r vertices that contribute there are suppressed by an additional m n /M [121 110j . 
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In general, to be in line with the Goldstone theorem, it is a necessary requirement that counter terms 
either contain a derivative acting on the pions, or scale with m^.. The latter is a consequence of the 
well known relation oc rh, where rh denotes the average mass of the light quarks. To be consistent 
with chiral symmetry the chiral Lagrangian is only allowed to have terms analytic in rh. To absorb the 
sensitivity to the wave function discussed in the introduction would require a counter term that neither 
contains powers of m, nor a derivative acting on the pion field at variance with chiral symmetry. 
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